Measure of genuine multipartite entanglement with computable lower bounds 
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We introduce an intuitive measure of genuine multipartite entanglement which is based on the 
well-known concurrence. We show how lower bounds on this measure can be derived that also 
meet important characteristics of an entanglement measure. These lower bounds are experimentally 
implementable in a feasible way enabling quantification of multipartite entanglement in a broad 
variety of cases. 



PACS numbers: 03.67.Mn,03.65.Ud 



O 

(N 
C 

o 

(N 



43 ■ 

Or 
i 

+-> • 

c ■ 
ctf ■ 
=5 

cr 



> 

o 
o 

(N 



Introduction 



Entanglement is an essential component in quantum 
information and at the same time a central feature 
of quantum mechanics [HQ- Its potential applica- 
tions in quantum information processing vary from quan- 
tum cryptography Q and quantum teleportation Q to 
measurement-based quantum computing j5]. The use of 
entanglement as a resource not only bears the question 
of how it can be detected, but also how it can be quan- 
tified. For this purpose, several entanglement measures 
have been introduced, one of the most prominent of which 
is the concurrence [3, 0, @|. However, beyond bipartite 
qubit systems Q and highly symmetric bipartite audit 
states such as isotropic states and Werner states 0, Q 
there exists no analytic method to compute the concur- 
rence of arbitrary high-dimensional mixed states. For a 
bipartite pure state \ip) in a finite-dimensional Hilbert 
space rl\®rli = C dl g) C d2 th e concurrence is defined 
as @ C(|V>)) = v /2 ( 1 - Tr / 9 i) where p\ = Tr 2 p is the 
reduced density matrix of p = \%p)(ip\. For mixed states 
p the concurrence is generalized via the convex roof con- 



struction C{p) = inf 



{Pi,\->l>i)} Y,iPi c (\^i)) where the inn- 
mum is taken over all possible decompositions of p, i.e. 
p = This generalization is well-defined, 

however, as it involves a nontrivial optimization proce- 
dure it is not computable in general. The concurrence 
is a useful measure with respect to a broad variety of 
tasks in quantum information which exploit entangle- 
ment between two parties. However, considering mul- 
tipartite systems, a generalization of the concurrence is 
needed that strictly quantifies the amount of genuine 
multipartite entanglement - the type of entanglement 
that not only is the key resource of measurement-based 
quantum computing [10( and high-precision metrology 



11 [ but also plays a central role in biological systems 
12 > [l3| , quantum phase transitions [bj . fla ] and quan- 
tum spin chains jlfjj. Although many criteria detect- 
ing genuine multipartite entanglement have been intro- 
duced (see e.g. Refs. [I3412|)> there is still no computable 
measure quantifying the amount of genuine multipartite 
entanglement present in a system. There are only few 
quantities available for pure states (a set of possible mea- 
sures is given in Ref. |33() which, however, are in general 
incomputable for mixed states and corresponding com- 
putable lower bounds have not been found so far. In this 
paper, we define a generalized concurrence (analogously 
to a measure proposed for pure states in Ref. [34]]) for 
systems of arbitrarily many parties as an entanglement 
measure which distinguishes genuine multipartite entan- 
glement from partial entanglement. As a main result we 
show that strong lower bounds on this measure can be 
derived by exploiting close analytic relations between this 
concurrence and recently introduced detection criteria for 
genuine multipartite entanglement. 



B. Genuine multipartite entanglement 

An n-partite pure state \ip) £ H\ ® Hi ® • • • ® H r , 
is called biseparable if it can be written as \ip) = 
\ipA) ® IV'b), where \ip A ) e Ha = H jl ® ... ® Hj k 
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and \iPb) € Hb = Hj k+1 <g> . . . ® Hj n under any bi- 
partition of the Hilbert space, i.e. a particular or- 
der {ji,j 2 , ■ ■ ■ jk\jk+i, ■ ■■jn} of {1, 2, • • • , n} (for exam- 
ple, for a 4-partite state, { 1, 3 1 2, 4} is a partition of 
{1,2,3,4}). An n-partite mixed state p is biseparable 
if it can be written as a convex combination of bisepara- 
ble pure states p = ^2pi\y^i){ipi\, wherein the contained 

i 

can be biseparable with respect to different bipar- 
titions (thus, a mixed biseparable state does not need 
to be separable w.r.t. any particular bipartition of the 
Hilbert space). If an n-partite state is not biseparable 
then it is called genuinely n-partite entangled. 
If we denote the set of all biseparable states by 6>2 and the 
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set of all states by 1S1 we can illustrate the convex nested 
structure of multipartite entanglement (see Fig. [I}. 

A measure of genuine multipartite entanglement 




FIG. 1: Illustration of the convex nested structure of mul- 
tipartite entanglement. The set of biseparable states £2 is 
convexly embedded within the set Si of all states (£2 C <Si). 

(g.m.e.) E(p) should at least satisfy: 

Ml E(p) = OVp e <S>2 (zero for all biseparable states) 

M2 E(p) > OVp e Si (detecting all g.m.e. states) 

M3 E(£ i P ipi) < Y,iPiE(pi) (convex) 

M4 E(Alocc[p\) < E{p) (non-increasing under local 
operations and classical communication) 40] 

M5 E(U i oca i pUl ocal ) = E(p) (invariant under local uni- 
tary transformations) 

There are of course further possible conditions which are 
sometimes required (such as e.g. additivity), but this set 
of conditions constitutes the minimal requirement for any 
entanglement measure. For a more detailed analysis of 
such requirements consult e.g. Refs. [33l. |35| . 

C. Concurrence for genuine n-partite entanglement 

Let us now introduce a measure of multipartite entan- 
glement satisfying all necessary conditions (M1-M5) for 
being a multipartite entanglement measure. 
Definition 1. For n-partite pure states |^) £ Hi ®Hi S3 
■ ■ ■ <£> H n , where dimi^Hi) = di, i = 1, 2, • • • , n we define 
the gme-concurrence as 

C gme (M) := min J2(l - Tr(p\ )) , (1) 

7; £7 V T » 

where 7 = {ji} represents the set of all possible biparti- 
tions {Ai|i3i} of {1, 2, . . . , n}. The gme-concurrence can 
be generalized for mixed states p via a convex roof con- 
struction, i.e. 

C gm e(p) = r inf y2p i C gme (\^ l }) , (2) 

% 

where the infimum is taken over all possible decompo- 
sitions p = For example, for a tripartite 
pure state e Hi ® Hi <8> H3 there are three 
possible bipartitions 7 = {{1|2, 3}, {2|1, 3}, {3|1, 2}}. 



Consequently, the gme-concurrence is C gme (ip) = 
min{ v /2(l-Tr(^)), y/2(l - Tr(p<)), y/2(l - Tr(pi))}. 
The definition of C gme {p) directly implies C gme (p) = 
for all biseparable states (Ml) and C gme (p) > for all 
genuinely n-partite entangled states (M2). Convexity 
(M3) also follows directly from the fact that any mixture 
Xpi + (1 — X)p2 of two density matrices pi and p2 
is at least decomposable into states that attain the 
individual infima. As the concurrence of any subsystem 
was proven to be non-increasing under LOCC (see e.g. 
Ref. 0]), the minimum of all possible concurrences 
will of course still remain non-increasing, thus proving 
(M4) also holds. Furthermore Tr(p 2 ) is invariant under 
local unitary transformations for every reduced density 
matrix irrespective of the decomposition, which proves 
that also condition (M5) holds. For pure states the 
gme-concurrence is closely related to the entanglement 
of minimum bipartite entropy introduced for pure states 
in Ref. [34[. In contrast to the original definition using 
von Neumann entropies of reduced density matrices, we 
use linear entropies. In this way we can derive lower 
bounds even on the convex roof extension which had not 
been considered before. 



D. Lower bounds on the gme-concurrence 

As the computation of any proper entanglement mea- 
sure is in general an NP-hard problem (see Ref. HH), it 
is crucial for the quantification of entanglement that re- 
liable lower bounds can be derived. These lower bounds 
should be computationally simple and also experimen- 
tally (locally) implementable to be of any use in practical 
applications. Let us now derive lower bounds on C gme 
which meet these requirements. Consider inequality II 
from Ref. [l?} , which is satisfied by all biseparable states 
(such that its violation implies genuine multipartite en- 
tanglement) 

^Ip^il^,... ^jij - J2 ^/(<i>|n 7 ^ 2 n 7 |$) < o , 
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(3) 

where |<E>) is any state separable with respect to the two 
copy Hilbert spaces and II{ a } is the cyclic permutation 
operator acting on the twofold copy Hilbert space in the 
subsystems defined by {a}, i.e. exchanging the vectors of 
the subsystems {a} of the first copy with the vectors of 
the subsystems {a} of the second copy. A simple example 
would be II{i}|0i02) ® IV'1^2) = l^i 02) ® l^i^)- 
For sake of comprehensibility let us show how to derive 
lower bounds for three qubits and then generalize the 
result (in the appendix). If we consider a most general 
3-qubit pure state in the computational basis 

\ip) = a|000) + 6J0Q1) + c|010) + rf|011) 

+e|100) + /|101) + .<?|110) + , (4) 



the squared concurrences C 2 (p 1 ) = 2(1 — Tr(p 2 )) with 
respect to the three bipartitions read 

C 2 { Pl ) = 4\ah-de\ 2 +F 1: (5) 
C 2 (p 2 ) = A\ah-cf\ 2 + F 2 , (6) 
C 2 (p 3 )=A\ah-bg\ 2 + F 3 , (7) 

where Fi are non-negative functions. The following rela- 
tions thus hold 

C(pi) > 2\ah - de\ > 2\ah\ - 2\de\ , (8) 
C(p 2 ) > 2\ah - c/| > 2\ah\ - 2\cf\ , (9) 
C{p 3 ) > 2\ah - bg\ > 2\ah\ - 2\bg\ , (10) 

and finally 

min{C( Pl ),C(p 2 ),C(p 3 )} 
>2|a/i|-2max{|de|,|c/|,|6 5 |} 
>2\ah\ - 2{\de\ + |c/| + \bg\) =: B . (11) 

Now for any given mixed state the convex roof construc- 
tion is bounded by 

C gme (p)> mf Y^PiBi. (12) 

{Pi,\Vi/} 

l 

For the choice |$) = |000111) and the abbreviation 
Puvwxyz '■= (uvw\p\xyz), inequality © reads 

I[p, |000111)] = Ipoooml - VPooiooiPiioiio 
— VPoiooioPionoi — \/Piooioo Anion < 0. (13) 

Now 

2|pooom| < inf y^pi2\a,ihi\, (14) 

{Pi,lw} 

due to the triangle inequality and 

2VPooiooiPiioiio > inf }Pi2\bigi\ , (15) 

{p%M>i)} i 

holds due to the Cauchy-Schwarz inequality (of course 
for all the parts of the other bipartitions). 
This leads to a lower bound on the convex roof construc- 
tion 

C gm e(p) > 2I[p, 1000111)] . (16) 

As C gme {p) is invariant under local unitary transforma- 
tions, we can infer that indeed every 2I[p, |<I>)] constitutes 
a proper lower bound. By taking into account the set of 
all vectors {|$)} we can thus define a computable lower 
bound which itself has many favorable properties (satis- 
fying Ml, M3, M4 and M5): 

C gme {p) >max2il>,|$)]. (17) 
I*) 

As the lower bound is straightforwardly generalized (the 
structure of the proof essentially remains the same, see 
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FIG. 2: Contour plot of the lower bound max| $ ) 2I[p, |$)] 
on the gme-concurrence for the set of three-qubit-states p = 
ci\GHZ)(GHZ\ + c 2 \W}{W\ + 1 - c l- C2 l given by convex mix- 
tures of a GHZ state, W state and the maximally mixed state. 
The greyscale is related to the bound max| S ) 2I[p, |$}] varying 
from to 1 (where is white), while the blue region denotes 
states which are positive under partial transposition with re- 
spect to all bipartitions. The optimization over all {|$}} was 
realized using the composite parametrization of the unitary 
group (see Ref. [371]'). 

the appendix for details) , eq. (1171) is indeed a proper lower 
bound on (jTJ) for any n-partitc qudit state. 

Discussion 

The detection quality of our obtained bounds on the gme- 
concurrence is illustrated in Fig. [5] for the family p = 
Cl \GHZ){GHZ\ + c 2 \W)(W\ + i^p^l of three-qubit- 
states, where 

\GHZ) = -^(|000) + |111)) and 
v 2 

\W) = ^=(1001) + |010) + |100)) (18) 

are the well-known genuinely multipartite entangled 
GHZ- and VT-state, respectively. It can be seen that 
the bounds are non-zero for a considerable amount of 
multipartite entangled states, especially in the vicinity 
of the GHZ-state. 

In fact, our bounds are exact for GHZ-like states, i.e. 
states of the form \gGHZ) = a\Q')® n + /3|l')® n wherein 
|0') G Hi and |1') G Hi are arbitrary mutually orthogonal 
vectors. By expanding \gGHZ) in terms of |0') and |1') 
analogously to (0| one finds C(pA„ r ) — 2\af3\Vji, hence 
Gg m e(p) — 2\a(3\. In order to prove the exactness of 
the bound we choose \<j>) = |0')®"|l')® n for inequality II 
which then yields 21 [\gGHZ)(gGHZ\, |</»)] = 2ja/3|. In 
fact we already know from the results of Ref. [IT]]? that 
the inequality will detect a huge amount of genuinely 
multipartite entangled mixed states in arbitrary high 
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dimensional and multipartite systems. In all of these 
situations we thus also have a lower bound on the 
gme-concurrence . 

Experimental Implementation 

In order to be useful in practice, measures for mul- 
tipartite entanglement need to be experimentally 
implementable by means of local observables (since all 
particles of composite quantum systems may not be 
available for combined measurements) without resorting 
to a full quantum state tomography (since the latter 
requires a vast number of measurements, which is 
unfeasible in practice). The bound (fT?) satisfies these 
demands, as for fixed |$) its computation only requires 
at most the square root of the number of measurements 
needed for a full state tomography. Furthermore it can 
be implemented locally as explicitly shown in [l8j]. In 
an experimental situation where one aims at producing 
a certain state \ip), it is now possible to choose the 
corresponding \<p) and not only detect the state as 
being genuinely multipartite entangled, but also have 
a reliable statement about the amount of multipartite 
entanglement the state exhibits. Even if the produced 
state deviates from the desired states, the criteria are 
astonishingly noise robust (as e.g. analyzed in Ref. fl7jk 
as for example a GHZ state mixed with white noise is 
shown to be genuinely multipartite entangled with a 
white noise resistance of rj 57%. 

Conclusion 

We introduced a measure of genuine multipartite entan- 
glement, which can be lower bounded by means of one 
of the currently most powerful detection criteria. These 
bounds are experimentally implementable and computa- 
tionally very efficient, allowing to not only detect, but 
also to quantify genuine multipartite entanglement in an 
experimental scenario. This has grave implications on 
applications where genuine multipartite entanglement 
is a crucial resource (as e.g. in quantum computing 
[j| or cryptopgraphy [38j]) and might allow to give a 
good estimate of the relevance of genuine multipartite 
entanglement in other physical systems (as e.g . in 
biological systems [l3| or quantum spin chains [Taj). 
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Appendix 

Let us finish by proving the lower bound for the general 
n-qudit case. For the most general pure n-qudit state 
IV>) = En,i 2 ,(...), in Ci u i 2 ^...) t i n \iii 2 (- ■ ■ )i n ) the squared 



concurrences C 2 (p 7 ) = 2(1 — Tr(p 7 )) with respect 
arbitrary bipartitions (7) always take the form 

c2 (Pi) = 4|c o(...) Cii(...)i - c a ( 7 )C /3 ( 7 )| 2 + F 1 , (19) 

(20) 

where F 1 are non-negative functions (see e.g. Refs. [H, 
l39l | for details on how to calculate the linear entropies 
of arbitrary subsystems). For every bipartition 7 there 
exists one pair a (7) and ,8(7) that can be retrieved from 

{a(7),/3(7)}=7T 7 {00(---)0,ll(---)l} (21) 

where 7r 7 permutes every number from the subset defined 
by 7 from the first half of the joint set with the second. 
Thus 

C(p-y) > 2|c 00 (... )oCii(... )i - c a ( 7 )C0( 7 ) I (22) 

will hold also for every 7. Now for the GME-concurrence 
we can infer 

min{C(p 7 )} > 2|c 00 (...)oCii(...)i| - (V \c a{l) Cf] h) \) =: B. 
7 

7 

(23) 

Now for any given mixed state the convex roof construc- 
tion is bounded by 

C gme (p)> inf V ft B t . (24) 

{piM'i)} 

For the choice |$) = |0)®" ® |1)®", inequality © reads 

= |A)0(-)011(-)i| - ^2 V^«(7)«(7)P/3(7)/3(7) - °' 
7 

(25) 

Now 

2|poo(-)on(-)il < r inf ,,y]p« 2 l c oo(-)o c ii(-)il' 

(26) 

due to the triangle inequality and 

VM7M7)P/3(7)/3(7) ^ , in f, , 5Z Pi2 l C «(-r) C /3(7)l > ( 27 ) 

(28) 

due to the Cauchy-Schwarz inequality. 
This leads to a lower bound on the convex roof construc- 
tion 

C gme (p)>2I[p,\$)]. (29) 

And again due to the local unitary invariance of C gme (p) 
this proves our lower bound for all | $) . 
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